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Abstract 

Scattering and production amplitudes involving scalar resonances are known, according to Wat- 
son's theorem, to share the same phase d{s). We show that, at low energies, the production 
amplitude is fully determined by the combination of 5{s) with another phase uj{s), which describes 
intermediate two- meson propagation and is theoretically unambiguous. Our main result is a sim- 
ple and almost model independent expression, which generalizes the usual i^-matrix unitarization 
procedure and is suited to be used in analyses of production data involving scalar resonances. 
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In reactions such as tx^ti'ti^ or ti^K~tt~^ Dalitz plots of experimental 

data display a wealth of resonating states (see for instance jll, Q]), indicating that contribu- 
tions from final state interactions (FSIs) to observed quantities are important. In order to 
disentangle basic interactions from this kind of empirical information, a reliable theoretical 
description of the final three-meson system is required. As the full treatment of this problem 
is rather involved, one usually resorts to an approximation, which consists in assuming that 
one of the final mesons acts as a spectator. The limitations of this approximation are not 
well established and the issue is being debated . In this quasi two-body approach, theo- 
retical work becomes much simpler, since one deals with two-body FSIs directly related with 
elastic scattering. This idea underlies Watson's theorem ^, which was produced more than 
fifty years ago and states that the same phase 5{s) occurs in both scattering and production 
amplitudes. Important as it is, the theorem does not determine how scattering informa- 
tion is to be used in the interpretation of production data and one finds a multitude of 
mutually contradictory prescriptions in the literature. In particular, there are parametrized 
production expressions based on either (sin 5) or (cos 5) [71, Isl, y, [lOl] . In other cases, ex- 
pressions used in data analyses bear no connection with the scattering phase. For instance, 
the pioneering study of the cr(500) scalar resonance in the decay ti^ti^ti^ jlj was 

based on a trial amplitude of the form 

n 

where nr stands for a non-resonant background, the v4j^*^* represent the contributions of 
various resonances and the complex weights Cj are determined from fitting to the data. 
The main ingredient of the A^^^ amplitudes is a relativistic Breit-Wigner function, aimed 
at describing resonance propagation and decay. This expression makes no use of the tttt 
scattering amplitude. An obvious problem with this loose adoption of trial functions to 
interpret empirical data is that information about the position of the resonance pole in the 
complex energy plane becomes contaminated by unreliable assumptions. 

We tackle the scattering-production problem in the case of a scalar resonance. In order to 
simplify the notation, we refrain from writing angular momentum or isospin labels explicitly 
in amplitudes and phases. We assume that the scalar resonance R is coupled to a pair of 
pseudoscalar mesons PaPh which, in practice, represents systems such as tttt, tiK, KK, in 
the presence of (j(500), fi;(800) and /o(980) resonances. The meson masses are fia and fib. As 



we are mostly interested in relatively low-energy processes involving pions and kaons, in this 
work we follow both the conceptual and mathematical descriptions of resonances proposed 
in Ref. 



The amplitude T describing the elastic process PaPb PaPb must respect unitarity, and 
below the first inelastic threshold, it can always be written as ^ 

Tis)=S^ Hindis) e'^^'\ (2) 

where 5 is a statistical factor such that Pa = Pb ^ S = 2, Pa ^ Pb ^ S = 1, and 
p{s) = — (/ia + /if,)^/s][l — (/ia — Unitarity is implemented automatically in 

this parametrization, by means of the real function which encompasses the full dynam- 
ical content of the interaction. 

In decays of heavy mesons, such as Z) or 5, a part of the width may be due to the direct 
production of a resonance at the weak vertex. When this happens, FSIs become important 
and the interpolation between the decay vertex and the observed \PaPb) state is described 
by the subset of diagrams shown in Fig. [H involving both the hare resonance propagator 
and the unitarized elastic scattering amplitude 7|, |9|. This subset of Feynman diagrams is 
represented by the function n(s) and, for the sake of conciseness, referred to as production 
suhamplitude. 

At low energies, n(s) is designed to replace the Breit-Wigner function associated with 
the scalar resonance in Eq. ([1]) and has the advantage of exhibiting a clear relation to the 
scattering amplitude. We demonstrate, in the sequence, that it can be expressed in terms 
of the elastic phase 5{s) as 

tan(5(s) 



_ ^, , cos Sis) 

rorm 1 : n(s) = g 



m|j — s 



e''^'\ (3) 

tanu;(sj 

where g is the resonance-mesons coupling constant, is the nominal resonance mass, such 
that (5(s = m?pj) = 7r/2, and uj{s) is a meson loop phase given by 

, , 7rp(s) . . 

where, 

^ The amplitude T is relativistic and we employ the conventions of Rcfs. 0| and [l^ . 
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^L{s) = p{s) log 



^'"^Mog(^). (5) 



This is our main result. Simple manipulations allow it to be recast in two alternative 
forms, namely 



form 2: n(s) = , e'^(^), 
form 3 : U{s) = g ^^^^ e^^(^) 



(6) 



(7) 



where Ai{s) is a running mass and T{s) is a running width. It is important to note that 
these two functions are not independent, since both of them are unambiguously determined 
by the phase uj{s). 
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FIG. 1: Resonance (wavy line) propagation and decay into mesons Pa (dashed line) and Pb (dot- 
dashed line). T stands for the the unitary scattering amplitude, Eq. ([2]). 



We now demonstrate the results given by Eqs. ([3])-([7j). In the spirit of the Bethe-Salpeter 
equation, the dynamical content of a unitarized elastic scattering amplitude can be realized 
as a series involving kernels and propagators. This idea is represented in Fig. where the 
iteration of the kernel /C by means of a two-body mesonic Green's function yields the full 
amplitude T. 

In the case of low-energy meson-meson scattering, the problem is much simplified by 
the fact that the two-body propagator can be described by a renormalized^ analytic loop 
function, which we denote by VL{s). In order to construct it, we note that the finite part of 



The mesonic loop integral contains an infinite contribution, which is eliminated by means of a coutertcrm, 
chosen so that one has J = 7r/2 at a point s = sq. As sq can be determined by measurements of the 
phase shift, for convenience, we write Sg = fTi|j.. Further details can be found in Ref. 12 1, especially its 



appendix B. 
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the two-meson propagator is given by a function L{s) such that 

• for < S < (yUfe - /ia)^ 



L{s) = p(s)log 



a/ (/ia + /ife)^ - S + (/ia - flb^ 



for (/ift - /ia)^ < S < (/ib + /ia)^ 



S - a/ (/ia - /ifc)^ - S 



log 



/ia 



(8) 



Lis 



/^a + 



for S > (/ift + /ia)^ 

- (/ia - ^bY - \/s- {lJ.a + l^bY 



2 I S \Ha 



(9) 



L(s) = p(s)log 



a/s - (/ia - /ife)^ + a/s - (/ia + /ift)^ 



/^a 



(10) 



where A(s) = (s — — fi^Y — 4:filfil. The expression for L(s) is simphfied if fia = ^b and 
can be found in 12|]. Above threshold, the function Q{s) is written as 

1 



5167r2 

= R{s) + il{s) 



[L(s) - L(m2j)] +z53L(s)} 



(11) 



and, by construction, R{m?j^ = 0. The imaginary component is particularly simple and 
reads I{s) = —p{s)/{S 16n). The phase uj{s) entering Eq. (jl]) is defined as 



tantu(s) = I{s)/R{s) . 



(12) 



The representation of T given in Fig. 2 corresponds to a geometrical series'^ and yields 



T = /c + )C{-n)}c + /c(-fi)/c(-fi)/c + 



(13) 



K K + K K 



FIG. 2: Geometrical series for the elastic unitary PaPb scattering T-matrix. K stands for the 
kernel, Eq. (fT4 ~ 



^ The fact that this series can be summed was shown by Oiler & Oset in Ref . ^13] . It is a particular property 
of the interaction kernel involving pions and kaons. 
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In a theory containing a resonance as an explicit degree of freedom, it is convenient to 
factorize its s-channel denominator in the kernel /C and we write 



— s 



(14) 



where 7(5) is a function satisfying the condition 7(m|,) 7^ 0. It is important to stress that 
Eq. (fT4l) corresponds just to a definition, which is completely independent of dynamical 
assumptions. The specific form of 7(5) for vrvr scattering in the framework of a SU(2) 
implementation of chiral symmetry was discussed in Ref. 12[, but it may be ignored in the 
present work. Using Eq. (fT4ll into Eq. (fT3|) . one finds 



T = S 

p(s) 7W^(s) — s — z m/j r(s) 
where we have defined a running mass and a running width by 

M\s) = ml + ^\s) R{s) , 
7^(5) p(s) 



m/jr(s) 



S16TT 



and, by construction, Al^(m^) 



m|j. Comparing this result with Eq.Q, one has 
tand(s) 



(15) 

(16) 
(17) 

(18) 



On the other hand, the elimination of 7^(s) from Eqs. f|T6l) and f|T7|) makes the constraint 
between the running functions evident and allows one to write 

1 1 f_ tan5(s) 



m 



R 



1 



tanuls , 



(19) 



after using Eq. (fT2|) . This result plays an important role in our derivation, because it indicates 
that the running mass is completely determined just by the information contained in the 
empirical function 6{s) and the theoretically reliable phase uj{s). 

The production subamplitude n(s) can be evaluated directly from Fig. [1] and reads 



His) 



mo — s 



As stressed by Pennington 



this shows that production and scattering share a "universal 
denominator", namely (1 +0/0), which imposes the Watson phase to the production sub- 
amplitude and transmits the physical poles from one amplitude to the other. Using Eqs. 
([13]) and dH]), one finds Eq. ^ 



His) 



72(5) 



/C 



72(5) 



mRT{s) 



sin 5{s)e 



i5{s) 



(21) 



which is form 3. Using Eq. (fTS!) it can be transformed into form 2, Eq. ([6]). Finally, Eq. 
f|T9|) yields form 1, Eq ([3]). This concludes our demonstration. 

Another interesting relationship that can be derived from Eqs. (13]) and (1211) is 

IGvr tan (5 



1C = S- 



(22) 



p{s) (1 + tan 5/ tanu;) ' 
showing that the kernel is determined just by 5 and u. 

Results ([3])-([7]) extend previous ones available in the literature. Our unit arizat ion, based 
on the two-body Green's function, generalises the usual /T-matrix approach 0, [m]. In 
the latter, the particles propagating inside the loop are taken to be on-shell, which amounts 
to ignoring the real part of VL{s): making ^ — > in Eq. flTTl) . one has cotu; and 
M."^ ^ m\. The deviation of form 1 from the i^'-matrix result is quantified by the factor 
(1 + tan5/ tancj). 




0.3 0.4 0.5 0.6 0.7 0.8 
s^^^ [GeV] 

FIG. 3: Ratio of the result from the iC-matrix approximation to our Eq. ([3]), given by (1 + 
tan(5Q/ tanw)"^, for the vrvr system. 



In order to produce a feeling for the numerical importance of this discrepancy, we con- 
centrate on the TTTT system and adopt the parametrization for the scalar-isoscalar phase shift 
Sq{s) quoted by Colangelo, Gasser and Leutwyler (CGL) |l5|. We compute n(s)/(7 using Eq. 
([3]) and identifying m^j with their ^Jsq = 0.846 GeV, which corresponds to 5o(so) = 7r/2. 
We show in Fig. [3] the ratio of the result obtained in a i^'-matrix approach to ours and 
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note that the former is typically 20 — 30% smaller in the range < ^/s < 0.8 GeV. In 
Fig. Hlwe show results for the real and imaginary parts of Il{s)/g, as well as its modulus. 
The profile of |n(s)/(7| has a broad peak-like structure close to the cr(500) pole position^ 
quoted in namely, ^/sZ ~ 0.475 — z0.290 GeV. This agrees with both empirical peaks in 
production experiments [l|, 
Chiral Perturbation Theory 



■ Il7| and theoretical works dealing with pionic FSIs based on 
18, 191. 
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FIG. 4: Result for Il{s)/g obtained by means of Eq. ([3]) using 6q{s) from 



In summary, we have presented a prescription that allows the production subamplitude 
of a low-energy scalar resonance to be determined from the elastic scattering phase, supple- 
mented by another one, describing a two-meson Green's function. This prescription has a 
number of nice features, namely: 

1. It has very little model dependence. As far as dynamics is concerned, the only 
assumption made is that the kernel entering the Bethe-Salpeter structure shown in Fig. [2] 
can be considered as point-like at low-energies. This assumption, which is standard in the 
study of mesonic systems by means of effective theories, allows one to use the rather well 
known bubble-like two-meson propagator. 



The result for the a pole positfon of 15| is not far from the latest and more precise determination 16 1 
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2. It contains the ii'-matrix approach as a particular case. In that approach, corrections 
from the meson loops to the resonance mass are neglecteded. However, as shown in Fig. [3] 
for the vrvr system, the loss of accuracy induced by this simplification is of the order of 30% 
in the region of physical interest. 

3. It gives rise to a bump which is also visible in experiments for the irn case. 

4. It indicates, by means of alternative forms 2 and 3, that expressions involving both 
(sin 5) and (cos 5) are acceptable, provided the coherent running mass and width are used. 
In case other forms for these functions are adopted, consistency is lost and results may go 
astray. 

5. It shows that precise measurements of both scattering and production amplitudes can 
produce direct information about the meson-meson interaction kernel. The effects of the 
unitarization procedure, encoded in the Q function, can be separated from the low-energy 
dynamics of the scattering process, which is represented by the function 7^(s) in Eq. (IT^ . 
Using Eq. ( 122!) one learns that this function can be determined from 6 and u. This is 
important because 7^(s) shapes the denominator of Eq. (fT5l) and plays a crucial role in the 
extraction of resonance parameters by means of pole extrapolation. Last, but not least, Eq. 
f l22|) yields directly JC in terms of 6 and u. 

6. It complements Watson's theorem, by showing that scattering and production ampli- 
tudes share not only a phase, but also neat information about the meson-meson interaction 
kernel. 

7. It can be tested experimentally by simply replacing a Breit-Wigner function with n(s) 
inEq. ([I]). 

A more detailed analysis of production data involving vrvr and irK resonances using the 
formalism developed here is in progress. 
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